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Математическая модель вязкопластичной смазки подшипников скольжения 
с деформируемой опорной поверхностью 


К. С. Ахвердиев, И. В. Колесников, М. А. Мукутадзе, И. С. Семенко 
(Ростовский государственный университет путей сообщения) 


Даётся метод расчёта радиального подшипника с деформируемой опорной поверхностью на основе аналога 
уравнения Рейнольдса для вязко-пластичной смазки и уравнения Ламэ для случая «тонкого слоя». Предла- 
гается метод, позволяющий формировать точное автомодельное решение рассматриваемой задачи. Полу- 
ченные аналитические зависимости позволяют оценить влияние безразмерного параметра пластичности и 
упругогидродинамического параметра на основные рабочие характеристики подшипника. В результате уста- 
новлено, что с увеличением значений параметра пластичности и упругогидродинамического параметра, зна- 
чения несущей способности и силы трения возрастают. В предельном случае, когда упругогидродинамиче- 
ский параметр стремится к бесконечности, значение несущей способности подшипника и сила трения стре- 
мятся к соответствующим значениям для подшипника с жёсткой опорной поверхностью. В принятом в работе 
приближении значение несущей способности и силы трения прямо пропорционально параметру пластичности. 
Ключевые слова: опорная поверхность, подшипники скольжения, вязкопластичная смазка, деформация, 
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Введение. Как известно, применяемые в настоящее время жидкие смазочные материалы (масла) 
состоят из масляной основы (базового масла) и композиции присадок, придающих маслам необ- 
ходимый уровень функциональных свойств [1, 2]. Добавки полимеров с высоким молекулярным 
весом придают маслам вязкопластичные свойства. Работа подшипников, работающих на вязко- 
пластичных смазках, достаточно изучена [3—6]. Анализ этих работ показывает, что в выполнен- 
ных исследованиях опорная поверхность подшипников считается абсолютно жёсткой. 

В области подшипников с жидкостной плёнкой смазки появилось новое направление под- 
шипников с нежёсткой опорной поверхностью. Жёсткость такой поверхности имеет такой же или 
даже меньший порядок величины по сравнению с жёсткостью плёнки смазки. Подшипники с не- 
жёсткой поверхностью имеют явные преимущества по сравнению с подшипниками жёсткой опор- 
ной поверхностью. Эти преимущества — допустимость больших несоосностей и деформации ра- 
бочей поверхности, а также терпимость к присутствию посторонних частиц между рабочими по- 
верхностями. Кроме того, податливость эластомера под действием давления смазки приводит к 
появлению своего рода губы, ограничивающей утечку сдавливающей смазки из подшипника. 

Таким образом, в условиях уменьшенного смазкой питания подшипник с нежёсткой по- 
верхностью может сохранить гидродинамическую или гидростатическую плёнку смазки. В этом 
отношении он значительно превосходит подшипник с жёсткой опорной поверхностью. Целью 
данного исследования является оценка рабочих характеристик такого подшипника и, в последу- 
ющем, сравнение их с хорошо известными характеристиками подшипника такого же типа, но 
имеющего жёсткую опорную поверхность. 

Постановка задачи. Рассматривается установившееся движение вязкопластичной смазки в за- 
зоре радиального подшипника с податливой опорной поверхностью. Шип вращается с угловой 
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скоростью в, а подшипник неподвижен. В полярной системе координат (1; 9) с полюсом в центре 

шипа уравнения шипа и деформированного контура опорной поверхности запишутся в виде 
=, Г’ = +ес0$0 + аф(6). 

Здесь г — радиус шипа; и, — радиус подшипника; е — эксцентриситет; аф(0) — функция, харак- 

теризующая деформацию опорной поверхности подшипника. 


При оценке влияния деформации опорной поверхности подшипника на его основные ра- 
бочие характеристики ограничимся максимальным значением функции аф(6). Введём обозначе- 


ния ар(9)= ар"; при 0 =0,2л| а” = таха(). 
Основные уравнения и граничные условия. Будем исходить из безразмерных уравнений 
движения вязкопластичной смазки для случая «тонкого слоя» и уравнения неразрывности, кото- 
рые получаются из уравнений Генки-Ильюшина методом оценок 
0°и 4 би ди 
= Е + А, + = 
ог” @а9 0г 00 
Система уравнений (1) решается при следующих граничных условиях 
и=0, ъ=0 приг =1+150$0 





0 (1) 








2 
и=-1 5=0 приг=0; р(0)=Р(2п) = 2%. (2) 
р 
е . * . . ' ого 
Здесь п= 5; б=и +аф" -п; о, = бо, о, =®@ки — компоненты вектора скорости; р’= 52 
| 21,5" ь 
гидродинамическое давление, г’=/, +0; А= В безразмерный параметр, обуславливаю- 


1 
щий вязкопластичные свойства смазки; р, — давление питания; то — предельное напряжение 
сдвига. 
К уравнению (1) необходимо добавить безразмерную систему уравнений Ламэ для «тонко- 
го слоя» 
3-08 
дг*? дг* 
Здесь в области занятой упругим слоем размерные величины и’, и, связаны с безразмерными 





0, =0. (3) 


соотношениями 
И, = ИИ, И, = ИИ, = пб", 0, =Г,-П, (4) 
где и" — характерная величина компонента вектора перемещений; 2-п — толщина упругого 
слоя. 


В переменных (/ 9) и (г*, 09) уравнение недеформированного контура, прилегающего к 
смазочному слою, можно записать в виде 
гг . е е 
г= т 1с050 = А, (6), г" =п, со$6 = ^, (9), п= 5' По 5 (5) 
т 
Уравнение деформированного контура и внешнего контура упругого слоя прилегающего к жёст- 
кой поверхности подшипника соответственно запишется в виде 
г =1+16056 = ^, (9), г" =1 +1, со$0 = ^, (9). (6) 
Система уравнений (3) решается при следующих граничных условиях 
ди ди. 
0 = би н М г 


ог" нь) 9/15) дг* 


М 

















= -Р, в |. 
г” = (0) 


19 


=, (в) = м] 


Вестник ДУ, 2012. № 8 (69) 








где 
* *Я2 
о. ЧО ЧО ) 
рег.б, (1-а)иолуб, 








С. — модуль сдвига; а — постоянная Мусхелишвили; р =тахр, 9 е 0, 21] ; р— безразмер- 


ное гидродинамическое давление в смазочном слое радиального подшипника с жёсткой опор- 
ной поверхностью. 
Граничные условия (7) означают равенство касательных и нормальных направлений на 
недеформированной упругой поверхности подшипника, прилегающей к смазочному слою. 
Интегрируя первое уравнение смазки (3) с учётом граничных условий (7), будем иметь 
р... Р 
и, =--— Г" +-— И, (0 (8) 
г М М , ( ) 
Воспользуемся приближённой формулой 
|. (6) -^ (6) мч. 


Ги" =рь (0 


С учётом формул (5) и (6) получим 





р 
ем (9) 


: 
Из уравнения (9) с точностью до членов [+ для аф’получим следующее приближённое 


уравнение 


= (и -п). (10) 


Из формулы (10), как и ожидалось, следует, что отношение максимального значения деформации 
к радиальному зазору прямо пропорционально безразмерному максимальному давлению и обрат- 


но пропорционально упругогидродинамическому параметру М. При М > ® аф“ 0. 


* 


аф 


Точное автомодельное решение задачи (1)—(2) будем искать в виде 
оф оф : й 
и=—+(г,0 =-——+И (г,0), = (8), И =И(5), 
5+ (гв), 5 =-55 + (6), = Ф(Е), И = (Е) 


и-5 (8) (6), =. Ф. ен (11) 

Подставляя (11) в (1) и (2) придём к следующей системе обыкновенных дифференциальных урав- 
нений и граничных условий к ним 

"в, б"=&, б'-Е'=0 (12) 





+'(0) 0; $ (0)=0, 2(0)=- и()=0, 5(0)=0,5(1)=0; [9 (Е)@Е=0. (13) 


Решение задачи (12)—(13) легко находится непосредственным интегрированием. В результате 
будем иметь 








2 
где с, =-6, константа С, в дальнейшем определяется из условия Р(0)=А2п). 
1 2 


- (Е -=), 6 (5) &2 Е 1 а (14) 


Гидродинамическое давление с учётом (14) определяется из последнего уравнения систе- 
мы (11). С точностью до членов С(п?) для определения Р приходим к уравнению 
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аР 3 
—+А=-6|1- 215050 + -п* + 
р П 5 
Из условия периодичности гидродинамического давления в принятом нами приближении для кон- 


станты 6. получим следующее выражение 


2 
т 05 2. + 6, (1- 316050 + 31? + 31? с0529) (15) 


© =А+6-ЗАп? - 91? (16) 
Интегрируя (15) с учётом (16) будем иметь 
р =—6пп@ — ЗАП @ + 91? 9129 + А 5т20 + ^^ (17) 
р 


где р=р (9) при 9 =6,, где 6, является корнем уравнения 

—6 с0$ 6, — ЗАЙ со$6, +181? с0$ 26, + ЗАЙ? с0$26, =0, й=е/(^ -^). (18) 
Перейдём к определению безразмерной несущей способности и безразмерной силы трения. 
С учётом формул (14) и (17) для К, — безразмерной составляющей несущей способности и без- 


размерной силы трения получим следующие выражения 
Ю 2п 
=— = - | Ряп@а@ = 6пп + ЗАпп 


кр 0 


о и Ию а 
тия р я в. в 21+ 41 п - Ап+ ри 
й 





у 


(19) 





(6) `^, (6) 
Численный анализ аналитический выражений (19) с учётом (17) и (18) проводился при следую- 
щих значениях параметров 





А=0,6; к —г, =0, 04 мм; - = 0,02 =й; М > 10*, 103, 10%, 10°; 0<М <10°. 
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Рис. 1. Зависимость безразмерной несущей способности от Рис. 2. Зависимость безразмерной силы трения от параметра 
параметра М: 1 — аф” = 0, А = 0,6; 2 — аф' =0, А = 0,6 М: 1 —аф’ =0,А = 0,6; 2 —аф'=0,А = 0,6 


Выводы. Из полученных выражений (19) и зависимостей, приведённых на рис. 1 и 2, следует, что: 

1. Значение несущей способности и силы трения подшипника с податливой опорной по- 
верхностью меньше, чем у такого же подшипника с жёсткой опорной поверхностью. 

2. С увеличением значения упругогидродинамического параметра М, значение несущей 
способности и силы трения подшипника возрастает. При М —› ® несущая способность подшипни- 
ка и сила трения стремятся к соответствующему значению для случая подшипника с жёсткой 
опорной поверхностью. 

3. В принятом в работе приближении значение несущей способности и силы трения прямо 
пропорционально параметру пластичности А. 
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